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SUMMARY 


Experimental measurements of the attenuation of plane shock waves 
moving over rough walls have been made in a shock tube. Measurements of 
the houndary-layer characteristics, including thickness and velocity dis- 
tribution behind the shock, have also been made with the aid of new opti- 
cal techniques which provide direct information on the local boundary- 
layer conditions at the rough walls. Measurements of shock speed and 
shock pressure ratio are presented for both smooth-wall and rough-wall 
flow over lengths of machined-smooth and rough strips which lined all 
four walls of the shock tube. A simplified theory based on Von Karmen's 
expression for skin-friction coefficient for flow over rough walls, along 
with a wave-model concept and extensions to include time effects, is pre- 
sented. In this theory, the shock-tube flow is assumed to be one-dimensional 
at all times and the wave-model concept is used to relate the local boundary- 
layer growth to decreases in shock, strength. This concept assumes that 
local boundary-layer growths act as local mass-flow sinks, which give 
rise to expansion waves which, in turn, overtake the shock and lower its 
mass flow accordingly. 

The results show that while agreement of boundary-layer measurements 
is good for all shock strengths, and while agreement of shock -attenuation 
measurements with theory is good for all the shock strengths in the smooth- 
wall case, the agreement of attenuation measurements with theory for the 
stronger shocks is poor in the rough-wall case. Discussions are presented 
which qualitatively account for this discrepancy in terms of distortions 
of the assumed one -dimensional wave model. 

INTRODUCTION 


There are no existing theories or experimental data which could be 
used to determine what attenuation might be expected with shock waves 
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moving over very rough surfaces. Some previous investigations have shown, 
by means of schlieren photographs, distortions of the shock near rough 
walls as well as the growth of a rather thick boundary layer along the 
wall. Investigations have also been made of the diffraction of waves 
about various-shaped bodies. In general, these data are o nly qualita- 
tive with regard to the shock strength as a function of time, distance, 
and surface conditions. 

To represent this type of flow with an exact theory would, of course, 
be exceedingly difficult if not impossible. In order to construct any 
theoretical model of this flow the problem must be highly idealized, in 
which case its applicability is usually quite limited. 

To transform the system directly to any stationary coordinates, 
such as flat -plate theory or impulsively actuated plate theory, would 
be a questionable procedure since in neither case is time allowed to 
sweep the plate, as it is in the actual case. The transfer of intelli- 
gence from the boundary layer to the shock is a process of great time 
dependence also. Three-dimensionality of the surface and the resulting 
wave diffraction pattern, as well as heat-transfer effects, all lead to 
great complication of the problem. 

In order to construct an experimental and theoretical model, even 
though highly idealized, from which comparisons might be made, the 
boundary-layer mass-flow sink concept from Donaldson and Sullivan (ref. l) 
was used, along with Von Karman's work on skin-friction coefficients for 
rough surfaces (ref. 2), as a starting point. From this starting point, 
which was first suggested by Donaldson (ref. l), and with corrections 
introduced to account for time of intelligence transfer from boundary 
layer to shock, an investigation was begun using the shock-tube facilities 
of the Langley gas dynamics laboratory. During the course of the inves- 
tigation a method using bullets was developed from which the boundary- 
layer height and velocity distribution at the rough wall could be 
determined . 


SYMBOLS 


a acoustic velocity 

A cross-sectional area; roughness constant 

B roughness constant 

Cf local coefficient of friction ' 
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mean coefficient of friction 
perimeter of inside of shock tube 
roughness height 

distance from leading edge of wall to measuring station 

mass flow 

Mach number, u/a 

velocity-profile power 

pressure 

Reynolds number based on x 

shock velocity 

time 

local velocity 

flow distance from leading edge of wall 
coordinate of boundary-layer height 

angle between flow direction and a line tangent to bullet wave 


5* boundary-layer displacement thickness 

0 boundary-layer momentum thickness 

5 nominal boundary-layer thickness 

0 shock strength, P 2 /P 1 

Subscripts : 

o high-pressure side of diaphragm 


1 


undisturbed conditions existing in tube ahead of shock wave 
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2 flow quantities in rear of shock-tube shock wave 

n Mach number of relative flow normal to wave from bullet 

T bullet wave relative to stream 

w flow quantities at wall 

b bullet wave relative to wall 

max maximum 


THEORETICAL METHOD 


In reference 1, the assumption was made that the boundary layer 
behind the shock wave, because of its growing displacement thickness, 
acts as a mass-flow sink and reduces the strength of the shock wave to 
correspond to the lower mass flow by the mechanism of expansion waves 
originating at the boundary layer and catching the shock, thereby 
weakening it. The same mechanism is herein used in calculating the 
shock attenuation due to smooth and rough walls for the case of turbu- 
lent flow, except that the time of travel of the expansion waves is 
included. The prediction of boundary-layer growth for both smooth and 
rough walls was taken from Von Harman' s work on skin-friction coeffi- 
cients (ref. 2). The method herein used also assumes a one -dimensional 
wave model as in reference 1. Wave diffraction effects are also neglected 
in this work, since in the one-dimensional case for small area changes 
the loss in wave strength is insignificant in departing from and coming 
back to the original area. 

/ / 

Von Harman's work shows essentially that in the case of very rough 
walls where the roughness heights are large compared with the thickness 
of a calculated smooth-wall laminar boundary layer, the boundary layer is 
a function only of the ratio b/k. That is, the dependence of rough-wall 
boundary layers is primarily on the roughness height and not on the 
Reynolds number, as is the case for smooth walls. 

For the smooth-wall case and turbulent flow (ref. 2) 


S = OOSxCf 1 / 2 


( 1 ) 


= 1.7 + 4.15 logio(R x Cf) 
Cf 


( 2 ) 
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where x, for the steady case, is the distance from the leading edge of 
the wall to the point under consideration. The value taken for x in 
this unsteady case is the maximum distance from the point where flow was 
initiated by the shock to a point where disturbances from this flow could 
be felt by the shock at a given location. Thus, the x distance depends 
on the location of the measuring station and the strength of the shock 
for a given length of wall. (See fig. 1. ) 
for x for the unsteady case is 


x up /U 2 + a2 

L a 2 \ Si 


Analytically, the expression 



where L is the distance from the leading edge to the measuring station. 
Equation ( 3 ) is plotted in figure 2 as a function of shock strength. 

With an appropriate assumption for the velocity profile, the boundary- 
layer displacement thickness is 


5* = pS OO 


Values of (3 are taken from reference 3 . The mass-flow loss is calcula- 
ble from 


Am = P 2112 AA 


( 5 ) 


Am » A A^ppUpj (6) 

Combining equations ( 5 ) and (6) and making the result nondimens ional 
gives 


a /P 2 ^2\ = P2 U 2 AA 
\Pp a -|_ J P]_ap A 


(7) 


With the approximation 


AA. = D5* 


( 8 ) 
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equation ( 7 ) may be combined with equations 


(4) and (8), yielding 


A 



. AA p S c 

Pl a l A 


(9) 


From a plot of the theoretical mass flow as a function of shock strength 
(fig. 3), the defect in shock strength can be found directly. 

For the rough-wall case, the problem resolves itself into finding 
the boundary-layer thickness as a function of x. From the boundary- 
layer momentum equation for a flat plate, 



( 10 ) 


Letting 


6 

0 


= 7 


(ID 


and 


de = 0 
dS S 


and combining equations (10) and (ll) yields 


2 db 

Cf = - — 
TJ dx 


or 


db = — Cf dx 
2 0 


( 12 ) 
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From reference 2, for fully developed turbulent flow over very 
rough plates , 


c f = 



(13) 


where A and B are constants depending on the form of the roughness 
(see refs. 2 and 4). 

Combining equations (12) and (13), letting 


dx = Up dt 


(14) 


and changing the dependent variable to 5/k gives 




+ B log e 


»f ® ^ 6. dt 

k/ k 2k 9 


(15) 


Integrating equation (15) yields 


6 

k 



- 2B(A - B) + 2B(A - B) log e | + B 2 log e 2 &■ 

K. K. 


u 2 t 5 
2k 0 


(16) 


Equation (l 6) is plotted in figure 4 for two sets of constants. 

Using the same x/L relationship developed for the smooth-wall 
case, there is a particular value of &/k corresponding to each shock 
strength. With a suitable assumption for the velocity profile, the loss 
in shock strength may be found from equation (9) and figure 3* A step- 
by-step procedure in calculating the shock attenuation would alter the 
results by less than 2 percent for the range of shock strengths 
investigated. 

The coordinate system used for the rough plates is shown in fig- 
ure 5* The position of y = 0 was taken at k/3, as a smooth wall at 
that height would give the same effective cross-sectional area. 
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The method of calculation for the determination of the velocity 
profile and thickness of the boundary layer by the bullet technique will 
now be discussed. (See fig. 6.) 

The method depends upon the assumptions that the wave strength just 
outside the boundary layer is maintained throughout the boundary layer 
and that there is no transfer of heat throughout the boundary layer. 

The first assumption should be valid for very weak waves regardless of 
the local boundary-layer conditions, while the validity for finite waves 
would depend upon density and velocity gradients in the boundary layer 
but would probably be evidenced by reflections of the wave. The second 
assumption is necessary in the absence of any knowledge of the boundary- 
layer temperature distribution but should again be valid in the lower 
subsonic range of stream velocities. 

From the first assumption the normal Mach number of the flow rela- 
tive to the bullet wave is constant, and therefore 


M n = Mrp sin 0*2 


u + % 

a 


sin a 


(17) 


where 


m t = 


^2 + % 
a 2 


(18) 


From the second assumption, and for the case of air near room temperatures , 


5a w 2 = u 2 + 5a 2 = U2 2 + 5a 2 2 


(19) 


Solving equations (17) to (19) simultaneously for u yields 


u D 5 sin 2 a 


u 2 M n 2 + 5 sin^a 


1 + \ M n 2 l^2 Z + 5 a 2 2 Mh 2 + 3 sin 2 a ^ 
' b si 


5 sinSal U)- 2 


5 sin 2 a 


u 

u 2 


( 20 ) 
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With the bullet velocity and stream-flow quantities in the shock 
tube determinable , the solution is then completed with measurements of 
a in the stream and in the boundary layer. 


TEST METHODS AND EQUIPMENT 


The shock tube was a rectangular tube of 4- by 7 tt -inch cross section, 

20 feet long, which could be evacuated or pressurized to 100 pounds per 
square inch. The arrangement of the shock tube and auxiliary equipment 
is shown in figures 7 and 8. The roughness was located on all four walls 
of the particular shock-tube section used. The roughness consisted of 
rows of 4-5° right pyramids machined at 45° angles with the flow and with 
the light beam through the window section. The pyramids were machined 
in plastic sheets which lined the tube walls, with the leading edges of 
the rough sheets faired and bonded to the tube wall. The sheets were 
made easily removable to facilitate the use of various heights of rough- 
ness. The location of the measuring gages in a small interchangeable 
section allowed data to be taken at various stations along the tube. 

Pressure and velocity measurements were taken of the shock wave at 
5- and 13-foot distances from the diaphragm. The 5-foot station deter- 
mined the shock strength before entering the roughness, while the 13 -foot 
station determined the shock strength as it came out of the roughness. 
Piezoelectric gages mounted in the tube wall, along with electronic 
chronographs, were used for measuring the velocity, while capacitance 
gages were used to measure pressures. 

The sheets of roughness produced a change of 13 percent in the effec- 
tive cross-sectional area of the 8-foot test section. Three heights of 
roughness were tested, pyramid heights of l/4, l/8, and l/l6 inch, but the 
effective -area change was kept constant. In addition, a set of smooth 
wood walls of equal area blockage was installed in the 8-foot section 
so that the effects of this blockage could be resolved by comparison 
with the smooth-wall data. 

The optical investigations were carried out in a small interchangea- 
ble section with 7g- “ hy l6-inch windows, and this optical section along 

with the schlieren equipment could be located at various positions along 
the shock tube. For optical studies, a strip of roughness was placed on 
one wall of the window section and the roughness section was not used. 

Two methods were introduced to study the boundary layer, a reflected 
shock technique and a bullet technique. The use of the reflected shock 
from the end wall of the tube permitted rough visual estimation of the 
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boundary-layer thickness by presuming that the change in curvature of 
the shock coincided with the outer edge of the boundary layer. 


The use of a bullet fired upstream and synchronized with the shock 
allowed for evaluation of both thickness and velocity distribution in the 
boundary layer by studying the changes in slope of the wave from the 
bullet nose as it propagated through the boundary layer. By regulating 
the speed of the bullet (adjusting the amount of powder used) the wave 
propagated into the layer was maintained at a speed slightly supersonic 
with respect to the wall velocity of sound. By firing small bullets 
relatively far from the wall the wave was held as close to a sound wave 
as possible. If the wave was finite, however, the picture could still 
be interpreted if the strength of the wave could be assumed constant 
through the boundary layer. 


The tests were conducted with Pp = 1 atmosphere for shock strengths 


p 2 

=1.6 to 

2. 5 

and 

Pi 


pi 


of 

Pp 

— = 2.5 

to — = 

5 . 8 . 


Pi 

Pi 


range of shock strengths and Reynolds numbers to be studied. 


RESULTS AND DISCUSSION 
Attenuation Measurements 


The experimentally determined shock attenuations at a station 
13 feet from the diaphragm are plotted as functions of shock strength 
for lengths of 8 feet and 4 feet of roughness, for roughness heights up 

to y inch, in figures 9 and 10. These results are cross plotted in fig- 
4 

ures 11 and 12 to show the dependence of the attenuations on roughness 
height. Typical pressure records are presented in figure 13 of shocks 


after traveling through 8 feet of i - inch roughness. Because measure- 

4 

ments of shock speed are more accurate than pressure measurements of 
this type, the pressure measurements were used mainly to show the quality 
of the pressure field behind the shock wave, rather than for measurement 
of the shock strength. Attenuation measurements using the wood filler 
strip are not plotted here since they were essentially the same as for 
the smooth case. Also, the results are not presented as functions of 
Reynolds number since no dependence on this parameter was noted in the 



= 1.6 x lo 5 


measurements for the range of Reynolds numbers covered 
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to 11 x 105 for P2/P1 of 1.6 to 2.5; ^7 = 1-1 X 10 5 to 6.3 X 10 5 
for P2/P1 of 2.5 to 5*sV 


Boundary-Layer Measurements 

Photographs illustrating the use of the reflected shock and the 
bullet wave to show the boundary layer are presented in figures 14 and 
15, as well as a photograph taken of an open jet used to confirm the 
method. The reflected-shock pictures are of academic interest only, 
since the distortion of the wave at the wall is a function of past events 
as well as present because the speed of the upstream wave is subsonic 
with respect to the wall velocity of sound. Since, in the case of the 
bullet wave, the wave speed (that is, bullet speed) is held supersonic 
with respect to the wall velocity of sound, the distortion of the wave 
represents only that due to propagation through the boundary layer at 
that point. A typical set of measurements of the change of wave angle 
in propagating through the boundary layer is presented in figure l6. 

By using the previously outlined method of equation (20 ) the veloc- 
ity distributions were calculated, and a typical one is shown in fig- 
ure 17. Figure 18 presents the results of all the boundary-layer thick- 
ness measurements by this method as a function of distance from the 
leading edge of the roughness. The averages of the velocity -profile 
power are also presented. The scatter of points in figure 18 at 
x = 3 feet is probably due to the influence of the weak cylindrical 
waves in the flow on the bullet wave. These disturbances are diffrac- 
tions which originated from the initial travel of the shock-tube shock 
wave over the roughness (see fig. 19)- In some of the photographs for 
x = 1 foot these waves were present in the boundary layer and those 
photographs were therefore disregarded; however, at x = 3 feet the 
waves were so weak that they were not always discernible. 


Theoretical Values of Attenuation 

From figures 14- and 18 and equation (l6) it can be seen that the 
selection of one set of values for the constants A and B in equa- 
tion (l6) will not exactly satisfy the measured values of 5 without 
using negative values of B. Because of the scatter at x = 3 feet the 
values of the constants were chosen to agree best with the values at 
x = 1 foot as well as the common assumption that 5=0 at x = 0. 

The values chosen were A = 3 and B = 3 even though the values A = 5 
an d B = 1.8 would not result in any serious change in the theoretical 
predictions. A constant value of velocity -profile power of l/5 was 



12 


NACA KM SL53D13a 


selected, based on the measured averages of close to l/4 and the fact 
that tabulated values were more readily obtainable for the l/5-power 
case (ref. 3)* Again, no serious change of predicted attenuation would 
result. 

It was then possible to calculate, using the method previously out- 
lined, all the theoretical attenuations for both the smooth- and rough- 
wall cases. The calculations were made for the smooth-wall case using 
a value of L = 13 feet and for the rough-wall cases using values of 
L = 9 feet and L = 5 feet, these lengths being the pertinent ones in 
this method, that is, the distance from the leading edge of the wall 
under consideration to the point of shock-strength measurement. The 
theoretical rough-wall attenuations consist of the sum of the attenua- 
tion calculated for the particular rough-wall length and that calculated 
for the smooth-wall length. The reason for adding these values in the 
rough-wall case is that the measurements are made under conditions whereby 
contributions to the attenuation are made by both the smooth and rough 
walls, and either contribution calculated under the assumptions used in 
this method is unaltered by the presence of the other; that is, the 
value of 5 max for the 13-foot smooth-wall calculation occurs ahead of 
the roughness strip, while the value of 5 mav for the rough cases occurs 

at a point where no flow particles with a history over smooth walls are 
present (see fig. l). 

A further reason for using the sum is that by representing a boundary 
layer by selecting its thickest point at a particular time, similar 
lengthwise distributions of boundary layer are assumed when comparing 
cases. In the case of mixed smooth and rough flow, however, it would 
not be suitable to use one or the other, but rather some combination of 
both should be used. The sum is merely the simplest combination, since 
they both certainly contribute to the attenuation. The theoretical 
attenuations are plotted in figures 9 to 12 as dashed lines. 

From examination of the experimental and theoretical results it is 
surprising (in view of the crudeness of the theory) to note that the 
theoretical predictions of shock attenuation are quite good for the 
weaker shocks at all values of roughness height and length, but for the 
stronger shocks they become progressively higher as the roughness height 
and length are increased. The differences in the cases of the stronger 
shocks are very great and are not to be expected in view of the boundary- 
layer thickness measurements . (That is, it should be remembered that in 
this theoretical method the variation of shock strength is mainly a varia- 
tion of effective flow length (fig. 4).) 

The measured values of boundary-layer thickness at both values of 
roughness length fit in very nicely with the theoretical values. Also, 
the measurements at 1 foot of roughness and at two greatly different 
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shock strengths confirm the theoretical prediction of invariance with 
shock strength alone. 

The most important fact shown by these results is that while agree- 
ment of boundary-layer measurements with theory is good for all shock 
strengths, and while agreement of shock-attenuation measurements with 
theory is good for all the shock strengths in the smooth case, the' agree- 
ment of attenuation measurements for the stronger shocks in the rough- 
wall case is poor. 


Distortions of Wave Model Not Allowed Under 
One -Dimensional Flow 

It is quite easy to conjecture as to quantities influencing shock- 
tube flows which are not accounted for in this theoretical model. How- 
ever, piecing the experimental evidence into the picture will eliminate 
many of the possible factors. For instance, if such factors as heat- 
transfer effects, boundary layers of cold gas flowing from the high pres- 
sure chamber, poor flow past a broken diaphragm, area changes due to a 
growing boundary layer or due to roughness, shock diffraction over rough- 
ness, and the influence of pressure gradients on the boundary-layer 
development are included in the theory, they would be expected to increase 
the predicted attenuations, produce second-order decreases, or produce 
attenuations that would not qualitatively vary as herein observed with 
shock strength and roughness height. 

One factor which might influence the theoretical predictions and 
which is not easily estimated is the effect of flow time on the boundary- 
layer development. (See fig. 20.) Since the theory used herein is taken 
from steady-flow flat-plate theory by substituting an effective flow 
length for the product of flow velocity and effective flow time, the 
boundary-layer particles in the shock case have not had the flow time 
which would exist in either the steady or the impulsive case. However, 
in view of the boundary -layer measurements at the 1-foot station (fig. 18) 
for greatly different shock strengths, and hence flow times, and in view 
of the good agreement of the smooth-wall attenuations, this factor must 
also be considered secondary for this comparison. 

Comparison of the agreement between the boundary-layer measurements 
and theory with the disagreement between the attenuation measurements 
and theory, in fact, seems to show in a very significant fashion that 
the marked deviations of experiment from theory are caused by factors 
dependent upon shock strength alone, that is, factors other than varia- 
tion of effective flow length with shock strength. The boundary-layer 
conditions at the wall are in agreement with the theory but this fact 
is not made manifest in the shock strengths observed. 



NACA KM SL53D13a 


Ill- 


Examination of equations ( 9 ) and ( 16 ) shows that the factors in 
the theory for rough walls influencing the shock attenuation are 


Attenuation = f 



k, M 2 , N, 



since 


and 


3 = f(M2, n) 


n = f(Mg, n) 


Boundary-layer measurements have ruled out all these parameters as 
causes of the discrepancy. 

The most logical source of error would seem to rest in the assump- 
tion that the rate of transfer of intelligence from the boundary layer 
to the shock is u + a. This rate is a function of shock strength only 
and would not always apply under the assumptions of the attenuation 
theory used herein, particularly at the stronger shock and larger rough- 
ness values. That is, in cases where the boundary- layer thickness (not 
displacement thickness) is of the order of the significant shock- tube 
dimensions, transfer of intelligence at the assumed rate of u + a would 
no longer be valid, since this quantity would no longer be propagated at 
this assumed rate. The boundary layer does, in fact, almost fill the 
tube for the stronger shock and larger roughness values. 

Figure 21 shows the parameters which would distort this assumed 
one -dimensional picture. The speed of transfer at the wall would be 
a 2 w > which is always less than S-^. The speed of transfer in the stream 

u 2 + a 2 i s .> °f course, always greater than S]_, so that the thickness 

of the boundary layer would determine the defect of transfer rate for 
a given shock strength. It can be seen, then, that the possibility that 
the speed of transfer will approach the value of S^_ might be very real 

at the thickest boundary-layer values. Figure 9 indicates those values 

of shock strength for the — - and i. - inch roughness heights, for 

k o 

L = 9 feet, at which the boundary layer actually fills the shock tube 
according to the theory. Boundary-layer measurements in figure 18 con- 
firm this fact. 
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The pressure records shown in figure 13 at a low and a high value 
of shock strength show qualitatively the fact that in the case of the 
stronger shock a train of expansion waves trails the shock. The curve 
of y'/g shown in figure 21 is included to show the effect on transfer 
of intelligence due to the two-dimensional aspect of the origin of the 
waves from the mass sink at the wall. It is not felt that this is of 
major importance , however, in this shock tube. 

It seems clear, therefore, that the defect In rate of intelligence 
transfer due to thick boundary layers is the primary cause for the devia- 
tions of experimental from theoretical rough-wall shock attenuations in 
this shock tube. 
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Figure 1.- Distance-time diagram used for theoretical calculations, 
showing contributions of smooth and rough walls to the measured 
attenuations . 











Figure 3.- Mass-flow rate across shock as a function of shock strength. 
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Figure 7 .- Arrangement of shock tube and instrumentation. 
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Figure 11.- Comparison of experimental and theoretical shock attenuations 
for 8 feet of roughness as functions of roughness size. 
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Figure 12.- Comparison of experimental and theoretical shock attenuations 
for 4 feet of roughness as functions of roughness size. 
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Figure 13. -Typical pressure records of a weak and a strong shock after 
traveling through 8 feet of i - inch roughness to show the difference 
in quality of the pressure field "behind the shock. 
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Figure 1^.- Boundary- layer thickness by 
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U) 0 = 3.5. 


reflected-shock technique. 
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(a) Distortion of wave from service 
air jet located over bullet path. 
(No shock-tube flow.) 



(b) 0.5 feet of rough wall; 0 « 1.7. 



(c) 2.5 feet of rough wall; 0 - 1.6. < d > x - 2 feet of rou S h wall J 0 * 3.9. 
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Figure 15 .- The bullet technique for finding velocity profiles and 

boun dar y-layer thickness for various lengths of rough vail and shock 
strength. 















Figure l8.- Comparison of experimental ‘boundary-layer thickness with 
^ | = ||a 2 - 2B(A - B) + 2B(A - B)log e | + B 2 log e 2 |]. 


theory . 


NACA RM SL53D13a 



Showing leading edge 


(b) 0 = 1.5. 
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(d) 0 = 4.6. 




-layer growth for various shock 


■e£l(366lS WH VOVM 






Figure 20.- Comparison of time histories of "boundary-layer flow for 
equilibrium with wall for three cases of flow from x = 0 to x. 
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Figure 21.- Parameters influencing one-dimensional concept of intelligence 
transfer from boundary layer to shock. 
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